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jo.  oralizn>d  clescribin  function  is  derived  for  a  dissymetric 
relay  with  D»  C,  co. iponencs  in  both  the  input  and  output  paths.  Methods 
of  evaluating  tie  D.  C.  components  are  discussed.  The  describing  func¬ 
tion  is  applied  to  the  analysis  of  step  and  disturbance  inputs. 

The  describing  function  is  odified  to  allow  computation  of  the 
D.  C,  and  sinusoidal  components  of  the  error  signal  when  a  ramp  input 
is  applied.  The  applicability  of  the  method  in  various  situations  is 
discussed. 

The  validity  of  tne  describing  function  analysis  was  investigated 
by  comparing  calculated  responses  with  actual  system  response,  as  simu¬ 
lated  on  the  analog  computer. 
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j:;  i 

i::**TODucrio:. 

In  receuu  ./cars,  tue  aescribin.,  function  ethod  has  been  viidely 
used  in  the  analysis  of  servo median  is.  is  containing  non-linear  elements 
(1),  The  basic  describing  function  technique  consists  of  "linearising" 
the  non-linear  portion  of  the  syste.i  by  assuming  that: 

1 .  The  input  to  the  non-linear  device  is  a  pure  sine  wave  of 

known  amplitude. 

2.  The  higher  harmonics  in  the  output  of  the  non-linear  device 

may  be  neglected. 

A  representative  describing  function  is  then  calculated  by  com¬ 
puting  the  amplitude  gain  and  phase  shift  between  the  input  and  the 
fundamental  frequency  of  the  fourier  series  expansion  of  the  output 
wave*  These  techniques  are  thoroughly  discussed  in  standard  texts  (2). 

In  some  syste  s,  however,  the  presence  of  the  non-linearity 
may  cause  a  D.  C.  co  lponent  to  exist  in  either  the  input  or  output 
path,  or  both.  This  is  particularly  true  with  dissymmetric  non- 
linearities  (3).  A  similar  situation  can  be  caused  in  many  symmetric 
systems  by  introducing  a  ra  ip  input.  In  either  case,  the  D.  C.  com¬ 
ponent  may  cause  a  change  in  the  magnitude  and/or  phase  of  the  funda- 
lental  frequency  of  the  output  from  the  non-linear  element,  and  it 
.may  result  in  a  significant  change  in  the  stability  characteristics 
of  the  syste.i. 

It  is  therefore  advantageous  to  define  a  generalised  describing 
function,  based  on  the  following  assumptions: 


1 


1.  i.i_  .  v.  i  :  •!.- linear  ele  :eni.  is  a  pure  sine  leave 

superi  pose  v  "  pure  2.  C,  component, 

2.  i.ie  higher  har.  onics  in  the  output  of  the  non-linear 

element  :ay  be  ne glee ted.  In  this  case,  the  output  of 
the  non-linear  ele  tent  lay  be  reduced  to  a  pure  sine 
wave  suneri.  .nosed  on  a  pure  J.  G.  component,  where 
botn  ter  :s  .ay  be  evaluated  from  the  fourier  series 
expansion  of  the  output  wave. 

The  generalized  describing  function  still  consists  of  the 
amplitude  gain  and  phase  shift  between  the  sinusoidal  component 
of  the  input  and  one  f unda  cental  frequency/  component  of  the  output 
wave. 

The  purpose  of  this  report  is  to  develop  the  generalized  de¬ 
scribing  function  for  a  dissy. metric  relay  and  to  investigate  the 
characteristics  of  servo mechanisms  containing  such  elements,  fhe 
report  includes  a  discussion  of  ;ethods  of  determining  the  magnitudes 
of  the  D,  C.  components;  an  analysis  of  system  stability  for  step  or 
disturbance  inputs;  and  an  analysis  and  discussion  of  the  system 
response  to  step  and  ramp  inputs. 
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J.-ki  jA  ±x 

1.  A.io::  OF  D.  C.  SIG  AL  L3V3L3 

2.1  General  Considerations 

In  all  cases  investigated  in  this  report,  it  trill  be  assumed 
that  the  block  diagram  of  the  servomechanism  can  be  reduced  to  the 
for  :  shown  in  figure  2.1,  there  G^s)  and  G2(s)  represent  the  linear 
portions  of  a  type  1  servo  :echanism,  and  G^  is  the  linearized  describing 
function  which  represents  the  non-linear  element  in  the  system. 


Figure  2.1  Block  diagram  of  a  servomechanism  containing  a  dissymmetric 
non-linearity. 


In  addition,  it  will  be  assumed  that  the  servo  loop  acts  suf¬ 
ficiently  like  a  low  pass  filter  so  that  only  the  first  harmonic  and 
the  D.  C.  component  of  O(t)  are  transmitted  around  the  loop.  There¬ 
fore,  for  steady  state  operation  with  sinusoidal,  step,  ramp,  or  dis¬ 
turbance  inputs,  we  can  approximate  the  input  to  the  non-linearity 
with  equation  2. 1 ;  and  the  effective  portion  of  the  output  from  the 
non-linearity  can  be  approximated  by  either  equation  2.2  or  by  equa¬ 
tions  2.3  through  2.5. 
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I(t)  =  '.-I  +  Xsin  uit 


2.2  0(t)  =  -~  +  a.cosujt  +  b.sinwt 

2  ‘  1 

2.3  0(t)  =  ~  4-  Csin(tnt  +  <}>) 


2.IJ  S  =  +  b,2 


2,5  =  tan-^  ~ 

b1 


It  v.'ill  also  be  assumed  that  G^(c)  represents  a  compensator 
having  a  transfer  function  of  the  form  expressed  in  equation  2.6;  and 
that  G2(s)  represents  a  rr.otor-load  combination  having  a  transfer  func¬ 
tion  of  the  form  expressed  in  equation  2.7. 


Gi  (s) 

=  k. 

7T(s  + 

zi) 

8 

! 

tF(TT 

Pi) 

where 

zi 

0  and 

Pi  t  0 

G2(s) 

=  k2 

TT  (a  +  Zj) 

1  s77(s  +  p  .) 

where 

zj 

0  and 

M  0 

The  basic  Laplace  transform  equations  for  such  a  system  ares 

2.3  £  =  ©i  -  ©o 

2.9  i(t)  =G1(s)(ei.e0) 


2.10  00  =  G2(s)0(t) 
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2 , 2  D,  C.  Level.  x.  i  i-r  Ut'-.'U'  Path 

In  the  case  of  ste->,  dis turbance ,  ran,  or  sinusoidal  inputs, 
the  D.  C,  co  ioonent  of  the  output  fro  :  the  non-linear  element  depends 
only  on  the  linear  pain  of  GpC5)  and  the  applied  input  function,  and 
is  relatively  easy  to  evaluate. 

The  equation  for  the  Laplace  transform  0(t)  can  be  obtained 
by  using  only  the  transfer  functions  of  the  linear  portion  of  the 
system  as  in  equation  2,11,  Since  the  steady  state  sinusoidal  com¬ 
ponents  have  no  effect  on  the  D,  C,  components  in  the  linear  portion  of 
the  system,  the  principle  of  superposition  can  be  used  to  eliminate  the 
transforms  of  these  components  when  solving  for  the  D,  C.  components 
only.  The  final  value  theorem  can  then  be  applied  to  obtain  equation 
2.12,  where  0(t)^,  ©■  - ,  and  I(t)c,  represent  the  Laplace  transforms  of 
the  non- sinusoidal  components  of  0(t),  ,  and  I(t),  respectively. 


2.11 


0(t) 


_  _  &o  _ 

u^IsT 


VCsT 


6i 


2,12 


s=o 


^2<s) 


I(t)d 
G-|  ( s  ) 


s=o 


In  order  to  simplify  the  following  calculations,  the  limits  of 
G^(s)  and  G^s)  as  s  approaches  zero  can  be  partially  evaluated  to 
get  equations  2.13  and  2.1':.  These  results  can  then  be  subs t it  its1 
into  equation  2,12  to  get  equation  2,15. 


2.13  G^s) 


s-*-o 


=  v1TT(s  *  Zi} 
77  (s  pi) 


s-*o 


_  ,  77  Zi  _  ^ 

‘TTVi  S 


5 


f  \ 

r.  *>  o 


77 


• J  '*  / 


-77 


S-vC 


- 

2~ 


s2 

r — 

“v2 


© 


id 


I(t),, 


V1 


s*»-  o 


for  a  step  input,  or  disturbance,  of  ,  agniuude  a,  the  La. place 
transfer  -!  of  is  given  oy  equation  2,1,  For  a  steady  state  analysis, 
the  initial  ti  ;e  is  in  .ate rial,  and  consequently  the  Laplace  transform 
of  I(t)d  can  be  represented  by  equation  2.1?  and  the  D.  C.  component  of 
C(t)  can  therefore  be  computed  by  equation  2.1". 


2.16  0.  =  Au(t) 

_  A 

1 

s 

2.1?  I(t)d=j 

2.13  22  =  ^! 

"  " 

=  0 

2  "v2 

s^-o 

a  similar  argument  can  be  applied  to  sinusoidal  inputs,  since 
in  this  case,  &id  is  zero. 

For  ramp  inputs  of  magnitude  3,  the  Laplace  transform  of  the 
input  is  given  by  equation  2.19,  and  the  D.  C.  component  of  0(t)  can 
be  evaluated  by  equation  2,20. 

2.19  =  3tu(t)  =  -2- 

s  * 
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For  an  n°n  order 


he  s  ler.  in  the  denominator  of 


Gp  ( s )  would  be  changed  to  s  1  Jquacion  2,21  can  be  derived  for  such 
a  system  by  following  a  line  of  reasoning  similar  to  that  used  for  the 
first  order  system,  Evaluation  of  equation  2,21  shows  that  the  d.c, 
component  of  0(t)  will  be  zero  for  both  step  and  ramp  inputs  for  any 
system  of  order  higher  than  one. 


ao  _  sn  ft 

T  ~  It 


2,21 


2.3  D«  C.  Levels  in  the  Inout  Path 

In  the  preceding  discussion,  there  was  no  limitation  on  the 
nature  of  the  non-linearity  represented  by  G^,  since  the  D,  G,  com¬ 
ponent  in  the  output  path  is  independent  of  G^,  In  general,  however, 
the  D,  C,  level  in  the  input  path  will  not  only  depend  on  the  nature 
of  the  system  input  and  the  characteristics  of  the  linear  portion  of 
the  system,  but  will  also  depend  on  the  nature  of  the  non-linear  ele¬ 
ment  , 

In  many  cases  the  solution  for  the  D,  C,  component  in  the  input 
path  can  be  extremely  complicated,  and  nay  result  in  transcendental 
equations  which  can  best  be  solved  by  graphical  means  (3), 

It  should  be  noted,  however,  that  in  many  cases  it  is  unnecessary 
to  solve  explicitly  for  the  D,  C,  level  in  the  input  path  in  order  to 
obtain  the  describing  function  curve.  This  fact  is  particularly  well 
illustrated  in  the  analysis  of  the  relay  servo  with  a  ramp  input,  given 


in  section  4,2,  In  other  cases,  such  as  the  analyses  for  step  inputs 
given  sections  3.2  and  3.3i  the  solution  for  the  D.  C.  level  indicates 
that  for  frequency  responses  analysis,  the  D,  C.  level  in  the  input  path 


ca..  us  ±3*  in'  li’i  !»■'*.  •  r-i  »/  is  treated  as  a  cy 

:etrical  co  pontirc, 


c:Ar  j:i  in 

TIU  RELAY  3ER'/C 

3.1  Discussion  of  the  Describin'-  Function. 

Appendix  I  contains  the  derivation  of  the  generalized  describing 
function  for  a  dissy 'metric  relay  with  the  following  characteristics: 

1 .  l’he  pull-in  voltage  for  a  positive  input  is  not  equal  to 

the  pull-in  voltage  for  a  negative  input. 

2.  The  drop-out  voltage  for  a  positive  input  is  not  equal  to 

the  drop-out  voltage  for  a  negative  input. 

3.  The  hysteresis  for  a  positive  input  is  not  necessarily  equal 

to  the  hysteresis  for  a  negative  input. 

4.  The  positive  relay  output  is  not  equal  to  the  negative  re¬ 

lay  output. 

The  input  vs.  output  characteristics  of  such  a  relay  are  shown 
in  figure  1-1  of  Appendix  I. 

Since  most  practical  relays  have  positive  and  negative  outputs 
of  equal  magnitudes,  no  attempt  was  made  to  analyze  a  relay  with  un¬ 
equal  outputs.  Therefore,  in  the  remainder  of  this  report,  it  will  be 
assumed  that  the  relay  output  voltages  are  of  equal  magnitude. 

3.2  The  Describing;  Function  for  a  Relay  vrith  no  Hysteresis. 

The  input  vs.  output  characteristics  for  a  dissymmetric  relay 
with  no  hysteresis  can  be  represented  by  figure  3.1. 

Equation  2.13  shows  that  for  a  servomechanism  of  the  type  shown 
in  figure  2.1,  the  D.  C.  component  of  the  relay  output  is  zero  for 
steady  state  operation  after  the  application  of  a  sinusoidal,  step,  or 
disturbance  input.  Combining  this  result  with  equation  1-15  of  Appendix 
I  gives: 


p 


Figure  3. 1  Transfer  characteristics  of  a  dissymmetric  relay  with  no 
hysteresis. 
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The  angle  relationships  of  equations  1-7  and  1-11  of  Appendix  I 
can  hot:  be  introduced  to  solve  for  the  D,  C.  co  roonent  of  the  input 
wave  and  to  solve  for  A  and  /?•->  as  shown  below. 

0. .  —  W  n  do  +  W 

3.3  smftj  =  _!_ -  =  sxnR  =  — 


3.^ 


di_  d2 

2 


sin  4  =  £LJ~„j=.  =  sin^0 

or 


JO  Sliir^j  —  — 


For  a  relay  with  no  hysteresis  the  expression  for  a.,  as  derived 


in  Appendix  I,  is 


3.3  a^  =  0 


hie  expression  for  b.  is  given  by  equation  I-l6  of  Appendix  I. 
co  Ibining  this  equation  with  equation  3.2  gives 


3.7  t.  =  ^.(cos B9  -!-  COS  $  )  =  z±cos£ 

7T  1  77  1 


4Y 

nr 


/'  -  s-i2ai  =f  i1- 


d1  J  d2 


2X 


the  equation  for  the  generalized  describing  function  is; 

-1 


q  P  n  _  1  /„  2  ,  1  2 

d  “  !T  /  a1  ‘  b1 


lan  p. 

£l 


For  the  dissymmetric  relay  with  no  hysteresis,  the  describing 
function  is; 
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for 


> 


Oh 


O 


d2 


Table  I  is  a  list  of  the  values  of  the  generalized  describing 
function  for  a  relay  vjix.li  d^  +  d£  =  5.2,  and  !  =  1 .  A  graph  of  this 
describing  function  is  shown  in  figure  3.2. 


X 

n 

r  -1 

Gd 

Of1  (db.) 

2.6 

0 

OO 

o° 

3.0 

.212 

4.73 

13.5 

3.6? 

.257 

3.^9 

11.5 

4.0 

.242 

4.13 

12.25 

5.0 

.217 

4.60 

13.2 

6,0 

.191 

5.23 

14.3 

7.0 

.169 

5.92 
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8.0 

.150 

6.65 

16,4 

10.0  | 

.123 

S.13 

1 _ 

17.2 

Table  I  Values  of  Gj  for  a  relay  with  no  hysteresis, 
d1  +  d2  =  ^.2 

Equation  3.9  applies  only  for  cases  in  which  the  relay  operates 
on  both  sides  during  each  cycle.  It  is  not  necessary  to  consider  one¬ 
sided  operation  at  this  point,  however,  since  this  type  of  operation 
requires  that  the  3.  C.  level  in  the  relay  output  .lust  be  different  from 
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zero.  Io  has  -reviousl,,  ^ee..  shown  Chat  for  a  servo .icchanis  1  of  the 
type  described  in  Chapter  II,  the  stead  state  D,  C,  component  in  the 
relay  output  vast  be  aero  for  sinusoidal,  step,  or  disturbance  inputs. 
If  the  D.  C,  level  in  the  output  is  allowed  to  be  different  from 
zero,  as  in  the  case  of  a  ramp  input,  then  both  one-sided  and  two- 
sided  operation  is  possible,  The  subject  of  ramp  inputs  is  treated  in 
Chapter  IV, 

3.3  The  Describing  Function  for  a  Relay  with  hysteresis. 

The  input  vs,  output  characteristics  for  a  dissymmetric  relay 
with  hysteresis  can  be  represented  by  figure  3.3  •  In  this  case,  the 
combined  results  of  equation  2,13  from  Section  2,2,  and  equation  1-12 
from  Appendix  I  gives: 


3.10 


ao  -  2L(* 
”2  2  IT  2 


3.11  /£+<,=^2  + 


Combining  the  angle  relationships  given  by  equations  1-3  through 
1-11  with  equation  1-13  of  Appendix  I,  we  obtain  equation  3.12, 


3.12  &1  =  ~ 


d1  "  Pi  +  y2  "  P2 


A 


If  the  hysteresis  is  the  same  for  both  sides  of  the  relay, 
d.j  -  p.j  =  d0  -  p,  ,  and  the  angle  relationships  can  be  evaluated  as 
follows : 

3.13  sin  A  -  sin  c*'.  =  Hi  =  H.2  ’2 

1  1  X  X 

=  sinfi  -  sino^, 
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Figure  3. 


lot  *  «(, 


to  t  3  IT  - 


lot  «2fT-^j 


3  Transfer  characteristics  of  a  dissymmetric  relay  with 
hysteresis. 


.  1 sin^j  -  si:  ^ 


2sin:;(^  -  o^cos-K^  +  <j) 

2sin-K^2  "  ^2^COS^^2  +  ^2) 


By  combining  the  results  of  equations  3 •  1 1 »  3-13*  and  3-  1if",  it 
can  be  shown  that  ^  ^  and  8^  - /?9<  The  D.  C,  component  of  the 

input  can  then  be  evaluated  by  equation  3.l6»  and  the  angles  can  be 
evaluated  by  equations  3*17  and  3. IS. 


5.15  sin/9.  =  —  _ _ 

A 


a2 


=  sin  4 


3.16  ■  =  d1  -  d2  ?1  -  ?2 

2  2 


3.17  sin^  =  sino<2  = 


_  ?1  +  P2 


2X 


.  1  sin^j  =  sin  4,  = 


2X 


=  R„ 


=  R, 


The  general  expression  for  b^  is  given  by  equation  1-1 4  of 
Appendix:  I,  Evaluating  this  equation  by  inserting  the  results  from 
equations  3.17  and  3. IS  gives  equation  3.19.  The  describing  function 
equation,  3.20,  is  then  obtained  by  substituting  equations  3.12  and 
3.19  into  equation  3.8.  The  describing  function  in  this  instance  also 
reduces  to  exactly  the  same  form  as  the  describing  function  for  the  cor¬ 
responding  symmetric  relay. 


3.19  b.  =  — (cos^i  +  cosO 

1  77  1  1 


=  f?<  J 1  -  %2  +  71  -  -v2  > 


16 


4>  =  i 


,  -1 


1  -  -R 


1  -  R, 


for  X  >  *  ?2 

2 


Table  II  contains  the  values  of  the  generalized 
cion  for  a  relay  with  +  dg  =  6  and  p.j  =  10. 

describing  function  is  shown  iir  figure  3.^. 


describing  func- 
A  graph  of  this 
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6,46 
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13.2 

7 
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5.75 

16.7 

10.1 

n 

.1372 

7.28 

17.2 

8.33 

9 

.1265 

7.9 

17,9 

7.15 

10 

.1163 

^00 

13.7 

6.26 

11 

.103 

9.27 

19.3 

5.62 

12 

.100 

10.0 

20 

5.1 

15 

.03018 

12.22 

21.7 

3.97 

20 

.  0624 

16.0 

24.1 

2.92 

30 

.0422 

23,7 

27.4 

1.92 

40 

.0316 

31.6 

29.9 

1.44 

50 

,025-4 

39.35 

31.8 

1.15 

Table  II 

Values  of 

G^  for  a  relay  with 

hysteresis. 

d1  +  d2  = 

6  and  P-]  +  P2  =  10 
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sa.ie  Tor  booh  sides  of  the  relay, 


the  tra:is coi.  eutal  e  .urtiuhs  which  .ust  be  solved  in  order  to  determine 
the  D.  C,  co  ponent  of  I(t)  become  m  .oh  more  complex.  They  are  difficult 
to  solve  with  Lrigono  :etric  relationships  alone,  but  it  is  relatively 
easy  to  solve  then  by  graphical  .  ethods  using  the  construction  shown  in 
figure  3,6.  This  construction  can  be  developed  in  the  following  manner: 


1 ,  Selected  a  prooractor  with  radius  R, 

2,  Calculate  the  follot.ong  quantities  for  each  value  of  X  to  be 


investigated. 


3.  Construct  figure  3*5  using  the  calculated  dimensions. 

4,  love  the  protractor  along  the  vertical  axis  until 

as  shown  on  figure  3.6.  This  satisfies 


equation  3.11. 

The  construction  guarantees  that  equations  1-3  through  1-11  are 

satisfied  if  and  only  ^ ,  ft  ,  and  are  selected  as  shown  in  figure 

3.6,  and  the  D,  C,  component  of  I(t)  is  established  by  the  center  of  the 

protractor,  figure  3.6  therefore  represents  the  desired  solution  for  the 
.-»•  ■» 

angles,  and  h  = 


The  generalized  describing  function  can  then  be  evaluated  vising 
equations  3.3,  3.12,  and  I- 1b.  Table  III  gives  the  values  for  such  a 
describing  function  for  a  relay  with  p^  -  3>  p?  =  7»  6.-J  =  2  and  dg  =  4, 
xhis  describing  function  is  plotted  in  figure  3.7. 
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Figure,  3.5 


Figures  3.5  and  3.6  Graphical  construction  for  solution  of  equations 

3.11  and  1-8' through' 1-1 1 . 
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lote : 

Values  for  and 

P  for  larger 

values  of  X  are 

approximately 

the  same  as  corresponding  values  in  table  II, 


Table  III  Values  of  the  generalized  relay  describing  function  for 
P-j  =  3*  P2  =  7,  =  2,  and  d2  =  4. 

3.4  Cxmarlson  of  Sy  metric  and  Dissymmetric  Relays. 

Equations  3.3  and  3.20  indicate  that  the  generalized  describing 
functions  for  dissy  metric  relays  are  exactly  the  sane  as  the  corresponding 
describing  functions  for  syr  .etric  relays,  provided  the  dissymmetric  relays 
have  equal  hysteresis  on  each  side.  A  servomechanism  containing  such  a 
dissymmetric  relay  would  therefore  respond  to  a  sinusoidal,  step,  or 
disturbance  input  in  a  manner  quite  similar  to  a  servo  containing  a  sym¬ 
metric  relay. 

If  the  servo  response  results  in  a  limit  cycle  or  sinusoidal  out¬ 
put,  the  output  oscillations  for  the  symmetrical  relay  would  be  centered 
about  the  position  commanded  by  the  D.  C.  portion  of  while  the  output 
oscillations  for  the  dissy  metrical  relay  would  be  offset  slightly  to 
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create  the  .ncessar  u,  Zt  me".  ach  signal  at  the  input  of  the  relay, 
rhus  the  effect  of  hiss;*  e try  is  to  introduce  a  steady  state  D,  C.  error. 
If  the  diss/  metric  relay  does  not  have  equal  hysteresis  zones  on 
each  side,  the  generalized  describing  function  nay  be  quite  different 
fro::  the  describing  function  which  would  be  obtained  by  assuming  equal 
hysteresis  zones  and  using  equation  3.20.  A  comparison  of  figures  3.4 
and  3.7  indicates  that  these  differences  will  be  relatively  minor  in  most 
practical  cases,  and  that  adequate  results  can  be  obtained  by  using  equa¬ 
tion  3,20. 

3.5  Verification  of  Proposed  Describing  F-.mctior.s . 

In  order  to  verify  the  describing  functions  previously  derived,  a 
position  control  system  consisting  of  the  dissymmetric  relay  and  third 
order  linear  components  was  postulated.  Linear  transfer  functions  with 
typically  low  pass  characteristics  were  selected  for  this  system. 

Figure  3*3  is  a  schematic  of  the  control  system  together  with  the 
analog  computer  simulation.  This  simulation  provides  for  operation  of 
the  relay  under  varying  degrees  of  dissymmetry  and  hysteresis.  (4) 

Figure  3*9  is  a  hichols  Plot  of  the  linear  transfer  functions  with  an 
arbitrary  gain  factor  of  20, 

Oase  I  Dissymmetric  Relay  'Without  Hysteresis. 

The  relay  xias  operated  without  hysteresis  and  with  a  constant 
total  dead  zone  of  5.1  volts.  The  plus  and  minus  pull-in  voltages  were 
varied  on  a  series  of  runs  to  obtain  various  conditions  of  dissymmetry, 
3tep  inputs  of  10  volts  and  15  volts  were  applied  to  the  system.  The 
gain  for  the  linear  portion  of  the  system  was  adjusted  to  different 
values  so  as  to  place  the  linear  curve  in  the  desired  position  on  the 
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fable  IV  is  a  siuriary  of  tie  results  obtained,  Typical  brush  re¬ 
cordings  fox-  several  values  of  system  gain  and  input  magnitudes  are  dis- 
played  in  figures  3.10  through  3.12,  which  illustrate  the  predicted  steady 
state  d.c,  error  levels  and  sinusoidal  characteristics. 
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Frequency 

1 

+2, 6, -2. 5 

QO 

•  >o 

4.6 

2 

4.7 

2 

2 

+3. 1 ,-2,0 

4,6 

2 

4.7 

2 

3 

+2.4.-2.7 

132 

6.0 

O 

3.0 

2 

4 

+4.0.-1.1 

132 

7.9 

2 

0.0 

2 

5 

+2.0.-3.1 

132 

3.1 

2 

0,0 

2 

6 

+3.5.-1.6 

132 

8.0 

2 

0 

CO 

2 

7 

+2. 7, -2,4 

81 
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63 

-.one 
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lone 
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Table 

17.  Experimental 

Results  for 

Dissymmetric 

Relay  Tithout  hysteresis. 

Case  II  Dissymmetric  Relay  with  Rystei’esis. 

The  relay  was  overate*  1  with  the  following  characteristics:  pull 
in,  +  .93  volts  and  -  3.  volts;  drop  out,  +  4.0  volts  and  -  1.93  volts. 
The  sane  general  procedure  was  followed  as  in  Case  I.  A  sunraiy  of  results 
is  shown  in  table  V.  typical  rush  recordings  are  displayed  in  figures 
3.13  through  3.13. 

A  co  iparison  between  the  predicted  and  observed  values  shown  in 
table  7  indicates  the  validity  of  the  proposed  describing  function  for  a 
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relay  wi  cl  hysteresis.  It  n  ....  .  e  ,,o  oeC.  that  a  small  difference  between 
the  observed  and  predicted  values  is  L-o  be  expected  since  the  experimental 
relay  characteristics  wore  nod  exactly  matched  to  those  used  for  predic¬ 
tions.  file  predicted  values  mere  based  on  figures  3.4  where  p-j  +  pg  =  10 
and  d-j  +  dg  =  6.  the  characteristics  of  the  experimental  relay  were : 

Pi  +  p 2  -  9.2  and  d^  +  d9  =  5.95. 
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4. 1  General  Gonsidera  .dons . 

..lien  a  relay  servo  is  subjected  to  a  step  or  disturbance  input, 
the  servo  '..'ill  either  eo  a  oo  rest  while  the  error  remains  in. thin  the 
lriits  of  the  dead  zone,  or  it  will  enter  a  limit  cycle  which  operates 
each  side  of  the  relay  during  each  cycle,  binilarly,  when  a  sine  input 
is  applied,  the  steady  state  error  signal  will  either  vary  id  thin  the 
limits  of  the  dead  zone,  or  it  will  cause  the  relay  to  operate  on  both 
sides  during  each  cycle, 

’.ihen  subjected  co  a  ra  ip  input,  however,  the  error  wall  never  re- 
rain  in  the  dead  zone,  the  servo  output  .ay  core  to  rest,  but  the  in¬ 
creasing  ramp  signal  in  uhe  input  channel  '..dll  always  cause  the  relay 
to  actuate  again.  the ref ore,  the  servo  rust  enter  a  Unit  cycle  when- 
ever  a  ra:rp  input  is  applied,  unless  the  ramp  signal  is  so  largo  that 
the  relay  does  not  drop  out.  this  limit  cycle  nay  involve  actuation 
of  both  sides  of  the  relay,  or  it  ay  involve  only  the  repeated  operation 
of  one  side. 

The  case  in  which  the  relay  operates  on  only  one  side  is  the 
simplest,  and  will  be  discussed  in  lection  4.2.  hen  the  relay  operates 
on  both  sides,  the  solution  must  be  obtained  by  graphical  methods,  as 
explained  in  section  4,3. 

4. 2  Describin';  Function  for  One-sided  Operation. 

.hen  the  relay  operates  on  only  one  side,  the  input  vs.  output 
characteristics  can  be  jorurayed  by  figure  4.1,  The  fourier  expansion  of 
the  output  wave,  ‘which  is  derived  in  Appendix  I,  becomes 
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Figure  4.1 


Transfer  characteristics  of  relay  operating  on  only  one  side. 
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7 r 


7T-  (°(  /3) 


'  a-1  -  £  (si  l#  -  sino() 

!  7 r 


=  -  (cos/9  +  cos  oO 


.Siei.  the  servo  has  a  ramp  input,  ©.  =  btu(t),  the  D.  C.  component 
of  the  relap  output  nave  is  pi von  bp  equation  2,20.  fhis  equation  can 
be  co  ibined  T.rith  equation  4.1  to  pet  equation  4,4 


2.20  fft  =  3 

^  kv2 


4.4  41  =  ior- 0^-/8) 

-^V2  ^ 

The  equations  for  the  generalized  describing  function  can  now  be 
obtained  in  the  following  manner,  without  solving  explicitly  for  the  D. 
C.  component  of  the  relay  input. 


4.5  a.,2  +  of  = 


IT 
--1 2 


7T 


=  2 


=  2 


2 

’'(sin2/?  +  sin2o<  -  2sin$sina0 
(cos  £  +  cos  o(  +  2cosificoso() 

(2  +  2coso(cosB  -  2sin  o(  sin  B  ) 

2 


1  +  cos(<\  + 


1  -  cos 


1  -  cos 


(7 T-  <£  -yfl)J 

22  n  1 

•vy  I 

“Vo  J 


4.6  tl  =  52bi  -  sliLgC  -  ^  1(  /  .  ^ ) 

b^  cos/9  +  cos  2 

4.7  ^  =  tan  -1  1± 
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ais  .  ,  ii  xs  convenient  to  introduce  the  trigonometric 
identity  given  by  equation  4.8. 


4, 3  sin£~  sinc(  -  2sin  i(^-o()cos  +  o() 

2  2 


=  2  sin  p  cos 


i r 

2 


bit 


in. 


V2 


=  2  sin  <j>  sin 


YK- 


■v- 


4.9  sin  p  =  Sn£=_ij£L* 


2sin  -.w 71 
Xhv2 


4.10  sin  /  = 


di  -  ?i 

2Xsta  JgL 
"v2 


The  generalized  describing  function  for  the  ramp  input  during 
one-sided  ooeration  becomes 


's-11  0d‘5f* 


1  -  cos 


23  7 T 


YIC 


v2 


0 


The  relay  will  operate  on  only  one  side  as  long  as  X  is  less 
than  a  critical  value,  Xc,  defined  by  equation  4.12. 

b  12  r  -  Pi  +  p2 

4.12  Xc  - 


1  +  sin  oi. 


If  the  relay  has  no  hysteresis,  there  should  be  no  discontinuity 
in  the  describing  function  as  the  servo  shifts  from  one-sided  operation 
to  two-sided  operation.  On  the  other  hand,  if  the  relay  has  a  signifi¬ 
cant  amount  of  hysteresis,  there  will,  in  general,  be  a  discontinuity 
when  the  shift  occurs. 
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contains  no  .niniaum  li.nl  o 


The  ex.  r«ssion  lor  u  .e  magnitude  of  3^ 
for  the  amplitude  of  the  sinusoidal  input*  The  phase  angle  expression, 
however,  indicates  that  the  describing  function  is  valid  only  if 


L'ote  that  if  the  denoininator  of  this  expression  is  equal  to  its  maximum 
value,  2,  this  restriction  lerely  requires  that  the  peak  to  peak  oscil¬ 
lations  of  the  sine  wave  must  extend  from  one  side  of  the  hysteresis  zone 
to  the  other. 

For  relays  with  absolutely  no  hysteresis,  the  describing  function 
predicts  absolutely  no  phase  shift  due  to  the  relay.  However,  if  the 
relay  contains  at  least  an  infinitesimal  amount  of  hysteresis,  as  all 
real  relays  do,  there  will  be  some  small  value  of  X  for  which  the  magni¬ 
tude  of  G^  is  defined,  and  for  which  the  phase  shift  will  be  90  degrees. 
Consequently,  the  effect  of  hysteresis  must  be  considered  to  be  much 
more  significant  in  the  case  of  ramp  inputs  than  for  other  types  of  inputs. 

The  describing  function  defined  by  equation  4.11  can  be  plotted 
in  non-dimens ionalized  form  by  defining  the  following  quantities 


4.14 


4.15  sin  f 


23sin  _JT 
—-v2 


- I—  Cm(s)  m(s) 

1  "  d1 


'  *1*  O  o 

4.1?  D?,  -  _J - 1 


The  critical  equation  for  stability  becomes 
4.13  3-d:.f1(s)l2(s)  =  4dnl(s)  =  -1 


The  plot  of 


Jd> ;  !db 


hZl 


can  then  be  used  as  che  ra.np  describing 


function  for  one-sided  operation  of  any  relay  by  simply  adjusting  the 


gain  of  the  linear  plot  by  the  factor 


?1  -  a. 


.  A  graph  of  the  de¬ 


scribing  function  of  equation  4.14  is  plotted  in  figure  4.2. 

.'he  naximun  value  of  2  for  which  the  describing  function  will  be 
valid  is 

4  10  7.  =  ,::c  =  Pi  +  ?2 

,  y  c  ppnq-  Tp1-v-T1Tn~-i-s'i?r^T 

— -  ^hc 

1  +  sin 

The  describing  function  plot  can  also  contain  loci  of  equal 
values  of  dp  =  Z(1  +  sin  o{ ) .  These  loci  will  pass  through  points  which 
correspond  to  the  naxi  an  value  of  2  for  which  one-sided  operation  trill 
occur  for  a  relay  with  the  corresponding  ratio  of  dead  zone  to  hysteresis. 

In  order  to  ootain  these  loci,  it  is  rost  convenient  to  first  plot 


Z  vs  Z(1  +  sin  o()  for  various  ratios  of 


This  requires  solvin: 


‘v2 


for  the  angle  $  using  equation  4.15  then  solving  for  the  angle  o(  using 
equations  4.4  and  4.7.  The  values  of  Z  can  then  be  obtained  for  constant 
values  of  and  these  points  cam  then  be  plotted  on  the  describing  function 
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Figure  4.5  Illustration  of  wave  shape  of  I(t)  in-  non-sinusoidal  region 
of  ranp  describing  function  plot. 
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curve,  v’Ow  .  _•  c  .Ic-lS  ol  3  vs  3(1  +  s±i<^)  for  various 

values  of  3  .  _  infer  afiou  o:.  these  curves  vras  used  to  obtain  the 

:Av2 

curve  for  =10,  x;hich  is  sno'.jn  on  figure  ’+.2. 

It  can  also  be  shown  that  there  is  a  liniwing  value  of  phase 


shift,  beyond  which  the  relay  input  can  no  longer  be  sinusoidal.  If 

3  >  ,  5YIC  ,  sinusoidal  operation  is  possible  only  if 
“2 


4.20  -0  5.  77- 

iAV2 

If  3  ^  .5IK  ,  sinusoidal  operation  is  possible  only  if 
v2 


4.21  -p  <  «j!L 

i^.r 

v2 

These  lira  ting  values  are  also  plotted  on  figure  4.2.  The  'mathematical 

proof  of  the  limitation  is  contained  in  Appendix  II. 

The  non-sinusoidal  characteristics  become  more  and  more  pronounced 

as  —  approaches  zero  or  unity,  and  as  p  approaches  90  degrees.  For 
mV2 

example,  for  3  =  .1  Ykv  and  p  =  81°,  the  solution  of  equations  4.4  and 

2 

4.7  requires  that  ck  -  1?2C  and  &  -  Q° ,  This  requires  a  wave  shape  of 
the  form  illustrated  in  figure  4.5.  The  corresponding  sinusoidal  wave  is 
shown  in  dotted  lines  for  comparison. 

4.3  Experimental  Procedures 

As  a  first  step  toward  the  determination  of  a  describing  function 

for  ramp  inputs,  a  series  of  investigative  runs  were  made  on  analog 

control  systems.  For  this  purpose  two  control  systems  were  simulated  with 

transfer  functions  ^  .  ,  and  ,  " _  ,  respectively,  hichols 

5Tb  +  4)  5T-  +  2)(d  +  4; 


Plots  for  these  systems  are  shown  in  figure  4.6. 
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*inaipsis  o.  w/.e  results  of  theso  preliminary  runs  together  x-rith 

the  theoretical  concepts  uhich  had  been  developed  up  to  this  point, 

indicated  the  significance  of  a  parameter  fomed  by  the  ratic  __JL  t 

Yh„ 


An  experiment  Xiias  then  designed  to  ’lore  completely  explore  this  feature, 

A  series  of  values  for  the  ratio  — —  ,  ranging  from  0.01  to  0.90, 

XaV2 


were  selected.  Several  runs  were  made  at  each  of  the  selected  values 
of  the  ratio,  by  varying  B  and  K  to  keep  the  ratio  constant.  The 
relay  x-xas  operated  at  a  constant  value  of  hysteresis  and  data  was 
recorded  for  both  the  second  and  third  order  systems.  The  reduced  data 
is  shox-rn  in  table  VI,  and  on  a  Nichols  Plot  in  figure  4.7. 

An  examination  of  the  data  plotted  in  figure  4.7  leads  to  the 
formulation  of  certain  conclusions.  Since  each  set  of  limit  cycle  condi¬ 
tions,  plotted  for  a  constant  value  of  %  s  seems  to  lie  on  a  single 

rKv2 

curve  it  appears  that  a  describing  function  curve  exists  for  each  value 
of  this  ratio.  The  fact  that  the  limit  cycle  conditions  obtained  for 
the  second  order  system  substantiate  the  pattern  established  by  the 
third  order  system  even  more  clearly  illustrates  the  existence  of  a 
describing  function.  The  nature  of  these  describing  function  curves  in¬ 
dicates  that  they  x-ri.ll  extend  across  all  phase  angles  from  -90°  to  -180° 
if  the  curves  are  extrapolated  beyond  the  region  of  plotted  data.  This 
agrees  x-iith  intuitive  reasoning  that  a  limit  cycle  must  exist  for  all 
type  I  systems  xrith  ramp  inputs,  unless  the  system  is  incapable  of  fol¬ 
lowing  the  ramp. 
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Phase 
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.5 

.01 

.924 

2.0 

103 

34.5 

.5 

.02 

.342 

4.45 

109 

36.2 

1.08 

1.47 

106 

26.7 

1.0 

.02 

1.515 

1.94 

111 

29.9 

.5 

.05 

.533 

2.40 

125 

23.5 

1.33 

0.95 

109 

16.5 

1.0 

.05 

.61 

3.67 

130 

28.2 

1.965 

1.21 

116 

19.0 

1.5 

.05 

.633 

4.8 

134- 

30.5 

2.40 

1.4 

121 

20.2 

2.0 

.05 

2.73 

1.56 

124 

21.3 

2.5 

.05 

3.14 

1.67 

128 

21.8 

.5 

.1 

.64 

1.65 

131 

15.7 

1.65 

.72 

112 

8.9 

1.0 

.1 

.796 

2.26 

140 

19.1 

2.325 

.86 

120 

10.9 

1.5 

.1 

.862 

2.95 

142 

21.7 

2.86 

.97 

125 

12.2 

2.0 

.1 

.932 

3.51 

146 

23.3 

d.p  1 

1.10 

129 

13.1 

2.5 

.1 

.96? 

4.70 

147 

24.6 

3.59 

1.29 

132 

14.1 

3.0 

.1 

1.03 

4.60 

149 

30.1 

3.92 

1.22 

134 

19.1 

1.5 

.3 

1.16 

1.47 

155 

8.3 

3.31 

.62 

129 

1.2 

3.0 

.3 

1.45 

2.10 

166 

10.5 

4.72 

.74 

140 

2.3 

1.5 

.5 

1.12 

1.08 

153 

4.0 

2.86 

.57 

125 

-1.8 

2.0 

.5 

1.21 

1.26 

156 

5.6 

3.40 

.55 

130 

-1.3 

2.5 

.5 

1.35 

1.43 

161 

5.8 

3.99 

.60 

135 

-1.1 

2.0 

.7 

.932 

1.11 

147 

5.7 

2.415 

.57 

121 

-  .1 

2.0 

.9 

.472 

1.16 

122 

11.7 

.967 

,62 

103 

6.6 

Table  VI.  Lata  obtained  in  experimental  determination  of  describing 
function  for  ranp  input.  See  figure  4.7  for  plot. 
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coincides  with  that  for  0,6,  the  curve  for  0,3  coincides  with  that  for 
0.7  etc, 

4.4  Verification  of  the  Proposed  Describing  Function 

Although  the  experimental  results  presented  in  section  4.3  con¬ 
firmed  the  existence  and  to  some  degree  the  nature  of  the  describing 
function,  another  experiment  was  designed  to  explicitly  verify  the 
theoretically  derived  describing  function.  In  addition  to  the  second 
and  third  order  control  systems  previously  simulated,  a,  fourth  order 


system  with  transfer  function 


K 


was  used.  In  con- 


S(S  +  1)(s  +  2) (S  +  4) 
ducting  these  runs  the  gain  of  each  linear  system  was  held  constant  so 

that  the  linear  curves  could  be  superimposed  on  the  plot  of  describing 

function  curves,  and  the  predicted  limit  cycles  were  obtained  from  the 

intersections.  The  value  of  the  ramp  input  was  varied  in  order  to  obtain 


different  values  of  the  ratio  JX 


YIC 


The  hysteresis  was  set  at  approxi- 


’v2 


mately  one  volt  and  held  constant  throughout.  Yichols  Plots  of  the 
three  linear  systems  are  shown  in  figure  4,6  and  the  nondimensionalized 
describing  function  curves  are  shown  in  figure  4.2.  The  predicted  and 
observed  limit  cycles  for  each  system  are  compared  in  table  711.  The 
Brush  recordings  are  displayed  in  figures  4.9  through  4.21. 

Since  some  of  the  predicted  limit  cycle  values  showed  dis¬ 
crepancies  of  25. j  or  more  from  those  observed,  it  is  felt  that  a  dis¬ 
cussion  of  error  sources  is  warranted  here.  First  of  all  the  inherent 
errors  of  any  describing  function  analysis  are  present,  that  is  the 
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-Uvy 

X 
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X 
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A 
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1 

.05 

.281 
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2 
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.562 
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,44 
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20.2 

3 

.20 

1.125 

2.89 

.61 
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.75 

25.8 

23.0 

4 

.30 

1.687 

2.79 

.748 

2.35 

.07 

15.7 

8.2 

5 

.50 

2.812 

2.73 

.350 

2.56 

.93 

7.6 

9.4 

6 

.70 

3.937 

2.75 

.722 

2.35 

.37 

14.5 

21.0 

7 

.80 

4.400 

2.85 

.572 

2.12 

.75 

22,0 

30.0 

Third 

Order 

System 

1 

.05 

.231 

1.74 

.379 

3.21 

.170 

84,0 

55.0 

2 

.10 

.562 

1.60 

.679 

1.77 

.56 

10,6 

17.5 

3 

.20 

1.125 

1 .49 

1.17 

1.32 

1.22 

11.4 

4.3 

4 

.30 

1.687 

1.40 

1.48 

1.32 

1.60 

8.6 

8.1 

5 

.50 

2.812 

1.36 

I.67 

1 .40 

1.35 

2.9 

10.4 

6 

.70 

3.937 

1.45 

1.39 

1.32 

1.60 

9.0 

14.3 

7 

.80 

4.400 

1.53 

1.07 

1.32 

1.22 

13.7 

14.0 

Second  Order  System 

1 

.05 

1.0 

1.21 

1.965 

2.52 

.75 

108 

62.0 

2 

.10 

2.0 

1.10 

3.31 

1.56 

2,4 

41.8 

27.5 

3 

.30 

1.5 

.62 

3.31 

.580 

3.37 

6,4 

1.8 

4 

.50 

1.5 

.57 

2.86 

.490 

3.00 

14.0 

4.9 

5 

.70 

2.0 

.57 

2.41 

.504 

2.40 

11.5 

0.4 

< 

O 

.90 

2.0 

.62 

.967 

1.25 

0.34 

101.5 

65.0 

Table 

VII. 

Recorded 

and  predicted  limit  cycle 

values  for 

■  ramp 

inputs. 
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describing  function  predicts  the  ratio  of  the  fundamental  sinusoidal 
components  of  the  input  and  output  signals.  In  many  of  the  runs  con¬ 
ducted  (large  and  small  ratios  of  )  the  input  to  the  non-linear 

iAv 

element  uas  decidely  non-sinusc idrl.  This  can  be  appreciated  by 
examining  the  Brush  recordings.  It  should  be  noted  that  the  value 
recorded  for  the  experimental  value  of  X  was  simply  ■§•  the  peak  to  peak 
amplitude  of  the  relay  input,  while  the  value  of  X  predicted  by  the 
describing  function  is  the  amplitude  of  the  fundamental  frequency 
comoonent.  Aopendix  II  shows  that  for  values  of  -JL-  near  0,1  or  0.9, 
and  for  values  of  fi  near  90°,  the  pull-in  angle  approaches  ISO0,  thus 
requiring  an  erctrenely  non-sinusoidal  input  signal  with  a  peak  to  peak 
amplitude  which  approaches  j  the  amplitude  of  the  fundamental  frequency 
component  as  illustrated  by  figure  4. 5. 

In  section  4.2  and  appendix  II  the  limits  for  sinusoidal  opera¬ 
tion  are  discussed  and  the  area  defining  these  limits  is  shown  on  the 
describing  function  curves  in  figure  4.2.  It  can  be  seen  from  the  table 
of  predicted  and  observed  values  that  the  errors  are  largest  outside 
the  area  of  sinusoidal  operation,  and  that  the  predicted  values  become 
increasingly  more  accurate  as  the  area  of  sinusoidal  operation  is  ap¬ 
proached. 

Another  consideration  in  evaluating  the  discrepancies  between  pre¬ 
dicted  and  observed  values  is  the  sensitivity  of  the  describing  function 
to  hysteresis.  The  analog  simulation  of  the  relay  employs  two  d.c. 
amplifiers  in  the  analog  computer,  set  to  high  gain  values.  It  was  not 
unco  mon  for  the  drift  in  these  amplifiers  to  cause  changes  in  the 
hysteresis  setting  of  several  tenths  of  a  volt.  In  order  to  gain  an 
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appreciation  Tor  :a  ,uxtuue  of  errors  introduced  by  s  .all  changes  in 
hysteresis,  a  representative  system  containing  fourth  order  linear  com. 
ponents  was  examined,.  The  results  of  this  examination  are  shoxm  in 
table  VIII  and  graphically  aispla  ed  in  figure  4.8,  The  extreme  sensi¬ 
tivity  of  the  amplitude  to  hysteresis  changes  is  evident.  It  is  also 
noteworthy  that  these  errors  are  essentially  invariant  with  the  ratio 


B 


In  conclusion,  after  careful  analysis  of  the  data  recorded,  it  is 
apparent  that  the  general  validity  of  the  proposed  describing  function  is 
proven.  Specifically,  when  intersection  with  a  describing  function  curve 
occurs  within  the  area  defined  for  sinusoidal  operation  the  amplitude 
and  frequency  of  the  limit  cycle  are  predictable  within  10m  error,  then 
the  intersection  occurs  outside  of  this  area  proportionately  greater 
errors  are  expected,  but  reasonable  and  useful,  estimates  can  be  obtained. 
4. 5  The  Ba.no  Describing  Function  for  Two-sided  Operation. 

If  the  sinusoidal  component  of  the  relay  input  is  greater  than 
the  critical  value  of  X  defined  by  equation  4.12  or  4.19*  the  system 
will  enter  a  limit  cycle  in  which  both  sides  of  the  relay  are  actuated 
during  each  cycle.  In  this  case,  the  pull-in  and  drop-out  angles  must 
satisfy  equation  4.20,  which  is  obtained  from  the  simultaneous  solution 
of  equations  2.20  and  1-12. 

4.22  +/S  -o(1  -/61  =  ZyJL 

1'<v2 

In  order  to  evaluate  the  describing  function  for  a  relay  with 

hysteresis  it  is  necessary  to  obtain  the  simultaneous  solution  of  equa¬ 
tion  4.22  and  equations  1-3  through  1-11,  This  solution  can  be  obtained 
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nError 

-j 

.3 

-30.8 

-4.6 

-16,3 

-1,2 

+15.7 

+1.2 

+33.6 

+2.3 

L 

•  • 

-30.2 

-2.2 

-16,3 

-1.1 

+15.9 

+1.1 

+35.1 

+3.3 

.5 

-31.3 

o  n 

-16,1 

-2.2 

+15.7 

+1.1 

+35.3 

+2.2 

.6 

-30.2 

-2,2 

-16.3 

-1.1 

+15.9 

+1,1 

+35.1 

+3.3 

.7 

-30.3 

-4,6 

-I6.3 

-1.2 

+15.7 

+1.2 

+33.6 

+2,3 

Table  VIII,  Errors  in  predicted  amplitude  and  frequency  resulting  fro:.: 

snail  relay  hysteresis  errors.  Zero  error  condition  is 
for  p-d  =  1,0,  Data  shown  is  for  a  system  with  fourth 
order  linear  components: 

2 

s(s  +'  V)'  "(s  +  2)  (s"  +  4) 
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Typies!  lir-xs'.i  req<prcllri£e  of  r&ure  responses. 
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ior  constant  values 


graphically  by  :lutiirg  o(  ~  0  o(  *,0  vse 

of  X,  The  value  c  vJ  car.  x  ■  1  be  So-.  insu  fcr  any  combination  of 

■q  rr 

and  X„  This  Infer. ration  can  an  turn  be  used  to  calculate  de- 

YK 

v? 

scribing  function  cur  ves  which  ari.ll  be  valid  for  taro-sided  operation, 
Unfoi’tunately,  this  nethoa  is  ncc  easily  adapted  to  a  non-dimension- 
alized  form.  Consequently  it  must  normally  be  calculated  individually 
for  each  relay  to  be  investigated. 

The  solution  js  slightly  less  complicated  for  a  relay  -without 
hysteresis.  The  transfer  characteristics  for  such  a  relay  are  shown  in 
figure  3. 1  „  In  this  case,  °(  -  /S  and  o(  =  A,,  Therefore  the  system 

of  simultaneous  equations  which  must  be  solved  in  order  to  evaluate  W 
consists  of  only  equation  4,23  and  equations  1-9  and  1-11, 


4,23  0.  -  /3X 


The  equations  for  the  coefficients  of  the  fundamental  components 
of  the  fourier  series  expansion  of  0(t)  reduce  to 


4,24  a  <5  —  C 


4,25  b  -  fl.  (  cus  m  cos/5-) 

•  rjj  .  £ 

The  solution  can  be  simplifies  further  by  noting  that 


4,26  cos  $  cos  0  v  =  2.  cos  f  (  B  <  +  £  ?)cosf-(^  -  0^) 


-  2  cosf(  |  9)cos 


2Yh, 


4,2?  b,  =  iicos  ;-(  $ ,  x  /30)cos 

rr 


2YK. 


v2 


?0 


Ail  eve  si  -lmr  oc  •  -ic  o...ri  je  ^  rained  for  a  system  containing 
an  ideal  relay.  In  this  case  tile  re  is  no  dead  zone,  and  ~  p 

Therefore  equation  4.27  reduces  to  equation  4.23  and  the  describing 
function  can  be  expressed,  by  equation  4.29. 


4.2 G 

4.29 


ud 


4Y 

rr  i 


COE 


It  should  be  noted  that  there  vail  be  a  discontinuity  between  the 
describing  function  for  one-sided  operation  and  the  describing  function 
for  two-sided  operation  whenever  there  is  a  significant  hysteresis  effect. 
Therefore,  if  the  system  is  operating  on  only  one  side,  and  a  small 
change  in  input  or  gain  causes  the  error  signal  to  change  so  that 
X  >  X„ ,  the  system  output  may  change  radically. 

It  is  possible  that  this  discontinuity  may  help  to  predict  some 
types  of  "jump  resonance"  found  in  relay  servos.  Suppose,  for  instance, 
that  the  system  is  subjected  tc  a  sinusoidal  input,  <9^  =  As in  ou^t,  such 
that  Auut  =  .9IXv2c  u, ,  Is  aueh  less  than  the  natural  frequencies 
of  the  system,  the  servo  will  behave  as  if  the.  input  signal  were  a 
slowly  changing  ramp  which  has  a  -laximuiri  slope  of  Au>0  at  u»0t  =  0  and 

decreases  to  zero  am  iot  ~  ^  . 

°  2 

how  suppose  the  system  characteristics  are  such  that  the  relay 
will  operate  on  only  one  side  when  ^0t  =  0.  Then  the  initial  operating 
point  will  be  determined  by  the  ramp  describing  function  for  one-sided 
operation.  The  operating  ooint  mil  tend  to  travel  along  the  linear 
curve,  with  the  instantaneous  operating  point  being  determined  by  setting 
B  =  Aui0cos  u>0t. 
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If  tho  linear  curve  crosses  the  critical  curve  of  for  the 
relay  used,  the  operating  point  will  eventually  reach  a  point  where  one¬ 
sided  operation  can  no  longer  occur,  and  a  two-sided  limit  cycle  will 
result.  During  the  next  portion  of  the  cycle,  the  servo  will  operate 
in  accordance  xjith  the  describing  function  for  two-sided  operation.  The 
servo  will  return  to  one-sided  operation  when  X  decreases  so  that 
Du  _ iL—  (1  +  sinc*p) ,  where  both  X  and  o<  „  are  determined  by  the 

Pn  * 't 

tx.’c-sided  describing  function. 

In  general,  the  magnitude  of  the  system  output  will  "jump"  to 
a  larger  value  when  the  relay  shifts  from  one-sided  operation  to  two- 
sided  operation,  and  will  "jump"  to  a  smaller  value  whenever  the  relay 
returns  to  one-sided  operation,  due  to  the  discontinuity  between  the  two 
describing  functions.  There  nay  also  be  a  change  in  the  frequency  of 
the  limit  cycle  when  the  shift  occurs.  Both  the  frequency  change  and 
magnitude  change  can  easily  be  predicted  once  the  describing  function 
curves  have  been  plotted. 

4,6  Sxtension  of  the  Rarro  Describing  Function  to  Other  Types  of  I~on- 
linearities. 

The  ramp  describing  function  can  be  extended  to  apply  to  other 
types  of  non-linearities.  In  general,  the  fourier  series  expansion  of 
the  output  wave  will  give  several  equations  of  the  following  types ; 

1 .  A  relationship  for  a 0  similar  to  equation  1-3. 

2.  Relationships  for  a*  and  b^  similar  to  equations  1-5  and  1-7. 

3.  Relationships  for  all  angles  involved,  similar  to  equations 
1-3  through  1-11. 
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jjLnce  e  .aAuiv-ji  c  ^  -jlies  for  any  non-linearity,  it  can  be 
combined  with  che  foarier  equivalent  of  a0  to  get  a  relation  between 
the  ramp  magnitude  and  the  angles  involved,  similar  to  equation  4.4, 

This  in  turn  can  be  combined  with  the  angle  relationships  to  solve  for 
the  d„c.  component,  .4  Once  b  is  evaluated,  it  is  a  relatively  simple 
matter  to  evaluate  the  equations  for  a-j  and  b-j  and  substitute  these 
values  into  equation  3=3  to  obtain  the  describing  function. 

It  should  be  noted,  however,  that  the  transcendental  equations 
involved  will  become  more  complex  as  the  nature  of  the  non-linearity 
becomes  more  complex,  and  the  solution  of  these  equations  may  become  ex¬ 
tremely  difficult.  In  addition,  other  types  of  non-linearities  may  have 
different  modes  of  operation  corresponding  to  the  one-sided  and  two- 
sided  modes  of  operation  for  the  relay  servomechanism,  and  there  may  be 
discontinuities  between  the  describing  functions  for  these  different 
nodes. 

Thus  a  complete  analysis  will  require  the  evaluation  of  a  separate 
group  of  describing  function  curves  for  each  node  of  operation,  along 
with  a  set  of  limit  curves  to  define  the  points  at  which  the  system  wall 
shift  from  one  mode  of  operation  to  another. 


73 


C.'-ir"  l'DA  V 

x-3co;4.*u;:dax'Io..o  ?oa  futuss  i.:v33tigatio::s 


During  the  course  ox  this  investigation,  several  interesting 
problems  were  noted  which  appeared  to  be  worthy  of  further  study.  These 
suggested  subjects  ares 

1 ,  A  comprehensive  study  of  the  effect  of  different  types  of 
dissymmetry  on  systems  containing  different  types  of  non- 
linearities  „ 

2.  A  more  exhaustive  analysis  of  the  ramp  describing  function 
for  a  relay  servomechanism  in  which  the  relay  operates  on 


both  sides.  This  analysis  should  include  an  investigation 
of  the  discontinuity  between  the  describing  function  for  one¬ 
sided  operation  and  the  describing  function  for  two-sided 
operation  to  determine  whether  or  not  these  describing 


functions  can  be  used  to  predict  certain  types  of  "jump 
resonance"  effects,  as  suggested  in  section  4,5. 

3.  A  thorough  study  of  the  extension  of  the  ramp  describing 

function  to  non-linear  devices  other  than  the  relay  would  also 
be  valuable  to  determine  the  effect  of  ramp  inputs  on  systems 
containing  such  elements.  An  outline  for  the  procedure  for 
such  an  investigation  is  given  in  section  4.6. 


74 


5I3LI0GRAPHI 


1 .  Kochenburger,  Rc  Je  Analysis  and  Synthesis  of  Contactor  Servo- 
aechanis-.s,  MTT  Servor.e char.is  .s  Laboratory  Report  ho,  1 6,  Seotenber 

1, 

2.  fhaler,  3.  J.  (  and  R„  3,  Broun ,  Analysis  and  Design  of  Feedback 
Control  Systems,  Me 3 rai  »  ill  Cook  Company,  Inc.,  Feu  York,  i960. 

3.  Langka.ner,  R.  P.3  and  t,  J,  biggins ,  Describing  Function  Analysis 
of  a  Type  1  Control  Syste.i  uith  an  Asy.r.ietrical  Deadzone  Yon- 
linearity ,  AIDS  Great  Lakes  District  Conference  Paper,  April,  1961. 

4.  Song,  Hui  3,  Simulation  techniques  for  "asic  Mon-linearities  in 
Servo; :echanis.r.s.  Master  of  Science  Thesis,  3,  S,  Maval  Postgraduate 
School,  Monterey,  California,  i960. 


75 


.  APPENDIX  I 


DERI  VAT  ICI'I  OF  THE  GENERALIZED  DESCRIBING  FUNCTION1 
FOR  A  DI3SYJ METRIC  RELAY 


Input  vse  output  characteristics o 


Figure  1=1 
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_ .  Fourier  series  expansion  of  output  wave 
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Tor  a  perfect  sir.e  in  rat,  the  angles  can  he  defined  as  follousi 


sin  = 


?1 


V/ 


~  '"O 


?1 


*  JH  -  '•/ 

1-9  sin  p4  -  — 1.— — -  =  R 
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Po  -  ".  i 

1-10  sin  qf  =  - - —  =  ’ 


a1 


1-11  sin  P2 


=  R, 
a2 


C.  The  describing  fun  -  cion. 


The  generalized  describing  function  is; 


If  both  sides  of  the  relay  are  actuated  on  each  cycle ,  equations 
1-4  through  1-7  can  be  used  in  the  evaluation  of  equation  3«S  without  re¬ 
strictions.  If  the  rale,  is  only  partially  actuated,  the  sane  equations 
nay  still  be  used  grovi-'c  •  certain  angles  are  defined  as  follows; 


If 


‘*o„ 


hi  .01 


n  =  1,2 
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Special  Cases  for 

"l  = 

Z2 

n  =  1,2 
n  =  1,2 
n  =  1,2 


if  Y-j  =  lr2  =  Y,  equations  1-3  >  1-5*  and  1-7  reduce  to 

:'12  ao  =W('/32  +  ^2  “  “  *1^ 

1-13  a1  =  ~(s in  /#  -|  -  sin  ^  +  sin  /6*  -  sin  ) 

1-14  =  ^(cos  /3  ■ j  +  cos  o(  ^  +  cos  /S  2  +  cos  o( 

If  Y-|  =  Y£  =  Y  and  the  relay  has  no  hysteresis,  then  =  By 
and  °( 2  =  A 2*  ana  equations  1-12,  1-13  and  1-1 4  reduce  to 

1-15  ao  =  fr  (  ^  2  “ 

3«e  a-j  —  0 

1-1 6  b-j  -  ~p(cos  A  7  +  cos  ^ ) 

If  Y.|  =  Y2  =  Y  and  the  relay  operates  only  on  the  positive  side, 
then  A  2  =  ^  o  ^  7<  The  subscripts  can  be  dropped  for  convenience  in 
writing,  so  that  equations  1-12,  1-13*  and  1-14  reduce  to 


4. 1  a0  =  |  (  TT  -  /S  -  *  ) 

4*2  a1  =  —(sin  /S  -  sin  o(  ) 

1  rr 
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=  ~(cos  0 


+  cos  o(  ) 
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Cr'iwt.’IC.:  0.  Ox.L-  07  J  FIDS 

.lien  one  relay  operates  on  only  one  side,  sinusoidal  operation 
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For  the  case  where  >  2T  ,  equation  II-2  reduces  to 
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